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^ ' Abstract: This talk is a review of news in low x physics which, I think, would 

ly^ ■ be useful for writing of the Monte Carlo codes. The following topics are discussed 

^ . here: (i) the next-to-leading order BFKLPomeron; (ii) two indications for shadowing 

Q>^ \ corrections (SC) in DIS from HERA data; (iii) matching of "soft" and "hard" 

^ \ photon - proton interactions; and (iv) survival probability for large rapidity gap 

^ ' ( LRG ) processes in hadron-hadron scattering and DIS. I hope, that our current 
understanding of these topics will allow us to narrow the gap between the MC codes 



and our microscopic theory - QCD. 
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1 Introduction 



Goals: The main goal of this talk is to share with you my understanding of the current 
situation in low x physics looking at this subject from the angle of possible improvement of 
existing Monte Carlo codes. Everyone knows that the typical MC code for deep inelastic 
scattering ( DIS ) contains three parts: the perturbative QCD cascade, the hadronization stage 
and the matching between "soft" anf "hard" processes. Since we have no solid theoretical 
understanding of the hadronization stage and the "soft" processes, there is a danger to write a 
MC code that will describe the experimental data, but in which the clear connection with our 
microscopic theory - QCD, will be broken. This is the reason that MC experts should follow 
all theoretical and experimental news in implementing in MC everything that we have learned 
both theoretically and experimentally. My goal is to provide them with a guide to the recent 
achievements in low x physics, which is a meeting point for all the complicated problems of 
perturbative and non-perturbative QCD. 

Topics: The choice of the topics, that I would like to cover here, is dictated by my per- 
sonal interests as well as by the considerable progress that has recently been made in their 
understanding. They are: 
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1. The next-to- leading order BFKL Pomeron; 

2. Two indications for shadowing corrections (SC) in DIS from HERA data; 

3. Matching of "soft" and "hard" photon - proton interactions; 

4. Survival probability for large rapidity gap ( LRG ) processes in hadron-hadron scattering 
and DIS. 

You can see that all above topics are closely related to our main problem- the matching between 
long distance ( non-perturbative ) physics and short distance ( perturbative ) approaches. 

Sources of information: This talk is a report on what I learned at: 

1. Theory Institute on Deep-Inelastic Diffraction - ANL, Sept. 14 - 16,1998; 

2. Workshop "Small x and Diffraction Physics" - Fermilab, Sept. 17 - 20,1998; 

3. 3-rd UK Phenomenology Workshop on HERA Physics - Durham, Sept. 21-25,1998; 

4. During animated discussions with J. Bartels, M.Braun, W. Buchmiiller, M. Ciafaloni, E. 
Gotsman, A. Kaidalov, Yu. Kovchegov, A. Kovner, J. Kwiecinski, L. Lipatov, U. Maor, 
A. Mueller, A. Martin, L. McLerran, D. Ross, M. Ryskin, G. Salam, M. Wiisthoff and 
many others, who shared with me their points of view on the topics that I am going to 
discuss; 

5. My own thinking on the subject, which has only partly been reflected in the published 
papers. 

2 The next-to-leading order BFKL Pomeron 



Why do we love the BFKL Pomeron? I think, we can list here three main reasons for 
our love to the BFKL Pomeron: 

1. The BFKL Pomeron |]l| is a high energy asymptotic in a perturbative QCD approach 
which we call leading log (1/x ) approximation ( LL(x)A ). In the LL(x)A we select the 
diagrams using as <^ 1, ds InQ^ ^ 1 but as ln(l/x) ~ 1. The LL(x)A is quite 
different from the leading log (Q^) approximation ( LL((5^)A , <^ 1, as ln(l/x) <C 
1, as 111 ~ 1 ) in which the DGLAP evolution equations were derived. We believe, 
that the BFKL Pomeron gives: 

(a) A guide for the matching of perturbative QCD (pQCD) with non-perturbative one 
(npQCD), at least for high enegy scattering amplitude in DIS; 
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Figure 1: The ZEUS'98 experimental data for forward jet production and their description 
in MEPJET MC, in CCFM evolution equation and in the BFKL approach ( BW 

(b) A possibility to get the "soft" Pomeron contribution by taking into account the 
npQCD corrections in the BFKL Pomeron. 

We have several examples which support our hopes 

2. The BFKL Pomeron is the only theoretical way to prove that at sufficiently high energy ( 
low X ) the density of partons becomes large. In other words, the BFKL Pomeron shows 
that the mathching of the perturbative QCD parton cascade with the non-perurbative 
high energy asymptotics, goes through a new non-perturbative QCD stage i.e. a high 
parton density QCD ( see, for example, review for details ). My personal opinion 
is that we do not need any theoretical arguments for a high parton density QCD, since 
the HERA experimental data show that the gluon density is high 0. However, the 
same experimental data we can be used as the experimental con&mation of the BFKL 
Pomeron; 

3. The new HERA data shows that the BFKL Pomeron is needed to describe the forward 
jet production ( see Fig.l). These data change the entire attitude to the BFKL Pomeron, 
which from a theoretical toy, becomes a tool for the description of the experimental data 
and, therefore, a part of DIS phenomenology which should be included in the MC codes. 
The question still remains could the experimental data of Fig.l be described without the 
BFKL Pomeron? You will answer this question better than me. However, Fig. 1 looks 
so impressive that it is difficult to believe that the agreement with the BFKL prediction 
is just a lucky coincidence. 
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Figure 2: 7*7* - interaction in the BFKL approach ( a ) and the time structure of the BFKL 
Pomeron ( h ). 

The map of disaster: It is well known that in the leading order the BFKL Pomeron leads 
to Regge-like asymptotics: 



atotiBFKL LO) oc e 



(1) 



where y — yo = ln(s/so) and r — Vo = In^Q'^/Ql ( see Fig.2-a for notations ). 
This simple asymptotic in the leading order is based on: 

1. The time structure of the parton cascade (see Fig. 2-b ); 

2. The separation between the longitudinal and transverse degrees of freedom; 

3. The fact that gluon emission leads to o" oc 

4. The diffusion in log of transverse momenta. 

Recently, the next-to-leading order corrections to the kernel of the BFKL equation has been 
calculated |Tl| and it turns out that these corrections are so large that they could change 
the entire understanding of the structure of the parton cascade in pQCD. Indeed, in the NLO: 



ujl < at g2 < 30GeV^ ^[^^ { see Fig.2-a ). It means that in a large kinematic 
region, the BFKL emission is suppressed, and only the Born term contributes; 

For Qq < 3 10^ GeV^ D < [|14| . It means that there is no diffusion in log of transverse 
momenta in the NLO BFKL. It also turns out that D < leads to oscillation of the total 



cross section which shows that the NLO BFKL has a serious pathology; 



In the NLO, the Regge-like asymptotics of Eq.( ^ ) is valid only in the limited range of 
energy y — yo = ln(s/ggo) < HH- For large values oi y — yo an extra term 

appears in exponent of Eq.( |l| ) which is proportional to {y — yo)^ ||l6l f^. Therefore, 



we cannot hope that in NLO BFKL we will get Regge-like asymptotics |jT6| [|T^ . 
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Even a slight glance at the difficulties, listed above, shows that the main properties of the 
LO BFKL are broken in the NLO. 



A possible way out - the main idea of resummation: The BFKL equation in the 
NLO can be written in the form: 



dG{y, q, go) 



S^'Kq-qo) + J d'q'K{q,q')G{y,(f,qo) , (2) 



dy 

where the NLO kernel K{q, (f) can be written as 

K{q,^) = as{q') { K^'^iq,^) + i^^^^(g',g') } = (3) 
^if,„„;(g-',g-') = ^{a5(go)K^^ + «|(go) if } . 



In Ref. JlTl it was found that the largest contributions to K^^^ with the Fadin-Lipatov 
choice of the energy variable y — yo = ln(s/ggo) is K^^'^ oc a| ln^(g^/g'^). One can see that 
such a term appears as a change of the energy scale in the normal double log contribution. 
Indeed, term a| ln(g^/gQ) ln^(s/ggo) leads to a| ln^(g^/g'^) ln(s/go) -contribution if we cange 
the energy scale from Sq = ggo to Sq = go- 

The above example shows that in the NLO we lost the very important property of the 
LO BFKL Pomeron, namely, matching the BFKL approach with the DGLAP one, for both 
g^ ^ Qo and q^ ^ g^. The main idea , suggested by G.Salam [|18| and by M.Ciafaloni |jT9|, 
is to replace the NLO kernel of the BFKL equation by the new kernel in which subleading 
- corrections to the BFKL equation have been resummed in a such way that 

TY-BFKL ^ T^DGLAP (*\ 

resummed ' \ / 

for g^ ^ gg and for q^ ^ g^. 

A possible way out - a toy model: Actually, the influence of the energy scale on the 
BFKL equation has been studied by Andersson,Gustafson and Samuelson pO| in their Linked 
Dipole Chain Model which reproduces the DGLAP double logs ln(l/x) ln(g^/gQ) for g^ ^ q^ 
and for q^ ^ g^. 



The key equation of this model 

dG{y,q,qo) agN^ f (Pi 



dy TT J nP 

with I = q — and 



{G(x',g'2) - G(x',g2)e(g-0 } (5) 



g^ , , g'^ 



y = ln(l/x) ; x = — ; x' = max(a;,a; • 

s q 

This equation is the BFKL one if x' = x. With a new definition of the energy variable x', the 
equation sums some of the next-to-leading order corrections to the BFKL kernel. 

Fig.3-a shows the Mellin image of the kernel for the BFKL case { x = x' , x^^^ )? for the 
next-to-leading order correction to the BFKL equation ( x^^~^^^) and for the kernel of Eq.( ^ 
) ( all orders). One can see, that the NLO gives a minimum ai u = which corresponds to 
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Figure 3-a Figure 3-b 

Figure 3: The Mellin image of the toy kernel (a) and the intercept of the toy-Pomeron (b ). 
See text for detail description. 



oscillation in the total cross section. However, the resummed kernel (all orders) has the same 
qualitatively behaviour as the LO BFKL kernel. In Fig. 3-b the BFKL intercept is plotted for 
the solution of Eq.( ^ ). One can see that the negative intercept for as > 0.2 in the NLO is 
irrelevant for the intercept of Eq.( |^ ), which is about twice smaller than the intercept of the 
LO BFKL equation, but is definately positive. 

Therefore, we can conclude from this toy model that the principle, formulated in the previous 
subsection ( see Eq.( ^ ), works and the resummation has a chance to heal our difficulties in 
the NLO BFKL equation. 

A possible way out - an example of resummation: 



G.Salam using the principle of Eq.( ^ ), suggested how to deal with unphysical logs 
of a| ln^(g^/gQ) - type. He demonstrated how to extend ( resum ) the NLO BFKL kernel so 
as to guarantee the concellation of unphysical logs to all orders and to provide the matching 
of the BFKL Pomeron with the DGLAP evolution equation. Unfortunately, the prescription 
for resummation is ambiguous and Salam suggested four realizations ( schemes) of his ideas. 
It should be mentioned that other schemes were advocated |]T2[ as well as the arguments for 



Salam's schemes 3 and 4 were suggested |[T^[^] [^p3[. In Refs. the criteria was 



suggested following the idea of Lipatov: the independence of the resummed result on the energy 
cutoff. This criteria is in the same spirit as the independence of the result on the factorization 
scale and it can be a good first step in understanding of the accuracy of the resummed BFKL 
Pomeron. 

Fig.4 shows the behaviour of the BFKL intercept (ul = ds x(|)) iii Eq.( )) and the width 
of the BFKL Pomeron ( D = asx"{^) )• One can see that the BFKL Pomeron intercept is 
positive in contrast with the NLO approximation. D is also positive, but only for schemes 3 
and 4. The second important conclusions that we can derive from Fig.4, is the similarity of the 
resummed NLO BFKL kernel with the toy model which has been discussed. 

Resume and recommendations: Let me share with you to-day ( 20.04.1999 ) my un- 
derstanding of the BFKL Pomeron standing. 
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Figure 4: The BFKL Pomeron intercept (ojl = ^'^'^ ^he width of the BFKL Pomeron 

(D = asx"i^,))- 

1. Large NLO corrections to the BFKL kernel left us without any theoretical predictions 
for both the BFKL Pomeron intercept ( ujl in Eq.( |T] ) ) and the BFKL Pomeron width 
( D in Eq.( |l] ) ). We can introduce both of them as phenomenological parameters in a 
MC code,, which should be extracted from the experimental data. The resummed BFKL 
kernel can be used to check how reasonable these extracted values of ul and D will be; 

2. For the running QCD coupling constant {as ) the Regge-like asymptotic of Eq.( || ) can 
be justified only in a limited range of energies. Therefore, it is better to use the numerical 
solution of the BFKL equation than the analytical expression of Eq.( [l| ); 

3. The main qualitative properties of the parton BFKL cascade remain the same as in LO 
BFKL, for the NLO resummed BFKL equations. This is the very important conclusion 
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for MC, which is the most striking theoretical result; 



It is better to use the Linked Dipole Chain ( LDC ) Monte Carlo [24] than the BFKL 
equation itself. However, we first need to check that the LDC MC describes the LDC 
model HQ. 



Two indications for shadowing corrections from HERA 
data 



HERA puzzle: The wide spread opinion is that HERA experimental data for > 1 GeV^ 
can be described using only the DGLAP evolution equations, without any other ingredients 
such as shadowing corrections, high twist contributions and so on ( see, for example reviews 
). On the other hand, the most important HERA discovery is the fact that the density of 
gluons ( gluon structure function ) becomes large in HERA kinematic region ( see Fig. 5 ). 

The value of the gluon density turns out to be so large that parameter 



Sn'^as xG(x,Q 
X - 



2Q2 



which determines the value of SC ||2^ ||2^ |7 



reaches unity 



(6) 



X ^ „ = 1 



2Ql{x) 



TT 



i?2 



(7) 



in HERA kinematic region (see Fig. 6 ) 

It means that in large kinematic region where k > 1 ( to the left from line k = 1 in Fig. 6 
), we expect that the SC should be large and important for the description of the experimental 
data. At first sight such expectations are in clear contradiction with the experimental data. 
Certainly, this fact gave rise to the suspicion ( or even mistrust ) that our theoretical approach 
to the SC is not consistent. However, the revision and re-analysis of the SC have been completed 
5U[ pT| [ P^ p5[ [^1^ with the result that k is the parameter which is reponsible for the 



value of SC. 

Therefore, we face a puzzhng question: where are SC in DIS? I am happy to tell you that 
we now where is the energy have two indications in the HERA experimental data that the SC 
are rather large and important. These indications are: 



1. xp - behaviour of the cross section of the diffractive dissociation { a£,£, ) in DIS; 

2. g2 . behaviour of Fa-slope ( ^^g^Q?'^ )• 

Here, we would like to discuss both phenomena and their relation to the SC. 

^ I am very thankful to L.Lonblad for discussions with me regarding the LDC MC. I realised that the MC 
experts have a problem in understanding why the LDC MC describes the experimental data worse than other 
approaches. I am very certain that as a result of such understanding will be a MC code which will include more 
from the BFKL cascade. 
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Figure 5: ZEUS ^B^j data for gluon structure function (xG{x, 



3.1 xp - dependence of ctdd 



Data: Both collaborations ( HI and ZEUS ) found that 

1 



X 



2 Af 



(8) 



where Ap 
• HI 



ap{0) — 1. The values of ap(0) are: 

ap(0) = 1.2003 ±0.020 (stat.)± 0.013 (sys.) ; 
ZEUSJ3§ : ap(0) = 1.127 ± 0.009 (stat.) ± 0.012 (sys. 



In Fig. 7 the ZEUS data on the Pomeron intercept are plotted together with the intercept 



the "soft" Pomeron It is clear that the Pomeron intercept for the difFractive processes in 
DIS is higher than the intercept of the "soft" Pomeron. 

Why is it surprising and interesting? To answer this question let us consider the 
simplest diffractive process - the diffractive production of the quark - anti quark pair with mass 
M ( see Fig.8 ) 
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Figure 6: Kinematic plot with lines k = 0.5, 1 and 1.5 . 
As it was shown ( see Refs. ||39l ) the cross sections of this process are proportional 



to: 



, , oc ^ X A ^ ^ / (9) 

dxpdt JQl 



for transverse polarized virtual photon and 

^^^.ML^ oc # X (10) 

dxpdt JQl kj 



for the longitudinal polarized photon. In these formulae we used the notation ( see Fig. 8 ): 
is the photon virtuality; M is the produced mass; xp = ^ ^^^^ , where s is the energy of 
virtual photon proton collisions; and (3 = g^rj^- 

The same factor also enters a more complicated process as, for example, a diffractive pro- 
duction of the quark - antiquark pair and one extra gluon (see Fig.8-b ). Indeed, the cross 
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Figure 7: ZEUS^d^ data on the Pomeron intercept for aj:,j:, in DIS. 




Figure 8-a Figure 8-b 



Figure 8: The simplest examples of the difjractive processes: ( a ) the production of quark - 
antiquark pair with mass M and (b) the production of quark - antiquark pair with one extra 
gluon. 



section of this process is proportional to |^ 



^/_^M!l^l±l±3. oc /"^ ^ X ln(MV4fcl) X x,G{x,,kl)f ^^^^ 

(XXpCit ^0 -L -L 

One can see that for the transverse polarized photon the diffractive cross section stems from 
the low values of the transverse momenta ( see Eq.( ^ ) and Eq.( ^ ) ) if we use for xG 



11 



the solution of the DGLAP equations. However, the situation changes crucially if we take into 
account the SC in expression for xG. For estimates of the value of the SC, we use the Glauber 
- Mueller approach for xG, namely 



xpG'^{xp,kl) = ^£ ^ /r ^/^^?{ 1 - e-'^"""""(-''^^)^(^*) } (12) 
vr xp X ^ r 

where S{bt) = e ^ . 

From Eq. ( |12] ) we can see that in the limit of low x 

xpG{xp,kl) =^ klR" ^\nK{x',-y) (13) 

Jx X k]_ 

if k^ <^ <5o(^)' where Ql{x) is plotted in Fig.6. Therefore, the substitution of xG^^ in Eq.( |] 
) and Eq.( |TT| ) leads to the situation that the integral over k^ is not infrared divergent, and 
the main contribution comes from the region of k^ of the order of Ql{xp). It means that 



daoD 

Xp 



Eq.( |l^ ) allows us to find the typical distances which are reponsible for the diffraction dis- 
sociation in DIS. Comparing the experimental values of X^jfiQ'^) with the ones calculated 
using Eq.( |12| ) (see Fig. 9 )0 , one can conclude that the typical k'j_ is not small, but rather 
k'j_ ^ 1 — 2GeV'^ 1^^,, if we take xp ^ 10"^, or even larger. 

It should be stressed that the diffraction dissociation for longitudinally polarized photon ( 
see Eq. ^ ) for quark - anriquark production comes from long distances and its cross section 
is proportional to 

d^t r ^ ^(^^' 4(1^) ) (^^) 

One of the striking feature of the experimental data, is the fact that the cross section of the 
diffractive dissociation in DIS and the total DIS cross section have the same energy dependence. 
Recently, an explanation of such a behaviour was suggested based on the SC contributions 142 



3| . Indeed, the total cross section and cross section of diffractive production are equal to: 

^totir?) = Jd'nJ^ dz\<iliQ^;n,z)\' atotirlx); (16) 
(^DDiYp) = Jd\J^ dz \^{Q^-rt,z)\^ {atotirlx) - (Tinirlx)} ■ (17) 

where atotirl^x) and o-DDijl^x) are total and inelastic cross section for dipole with distance 
between quark and antiquark, respectively. 

The solution of the unitarity constraint for dipole scattering gives: 



aUrlx) = 2Jd%,{l - e--. ^XM}; (ig) 



,{rlx) = Jd\{l - e-2-?^^^^^(M} ; (19) 



^Doing this comparison, we have to take into account that Xeff{Q^) = X{Q'^){Fig.9) — \ ^^"(i/x)^ • Qai^) 



in Fig.6 can be described by simple formula (5q(x) = IGeV'^ { ^)^-^'^ 
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where k^^^^^ 

After integration over z we get for transverse polarized photon 

1 



dz \^!T{Q'^-rt,z) 
Jo 



4 OL^Yn-^c 

3 7r2 Q2 



X 



^4 



(20) 



Using the expihct expression for S{ht) and assuming that the anomalous dimension of xG is 
not large 7 <^ 1 ( in this case we can replace xG{x, k^) by xG{x, Q^) in the integrals ) we can 
obtain El: 



atotirp) oc 2R^\n{K,) + C {1 + Kg) E,{k,) + 1 - e""' } ; 
cTDDirp) oc 2R^\n{Kg) + C (1 + Kg)E,{Kg) + 1 - e-^"} 
- { ln(2 + C ( 1 + 2 ) Ei(2 + 1 - e"2«. } ; 

where C is the Euler constant and Ei is the interral exponant [^. 



(21) 
(22) 



These formulae indicate that (Jtot ^ xG{x,Q^) and aon {xG{x,Q^)) for small <^ 1, 



but at ^ 1 we obtain that 



2i?2{ln(/s:,) + C}- 



(23) 



CTDoiYp) 



2R^\niKg) + C} - R^H2k,) + C} = i?Mn(«:,) = 1^^4.(24) 



This simple calculation shows that indeed at low x, where Kg is expected to be large, aoD 
and atot 0], should have the same energy dependence but accuracy of such a statement is 
controlled by ratio C/ln{Kg) -C 1. In Refs. 



and one can find estimates of how SC 
reveal themselves in real experimental situation. It turns out that they are able to describe the 
experimental energy behaviour of aon- 



3.2 The effect of screening on - dependence of the F2 - slope. 



Data: The experimental data ||4^ for the F2 - slope is shown in Fig. 10 ( Caldwell 

plot ). These data give rise to a hope that the matching between "hard" ( short distance ) and 
"soft" (long distance ) processes occurs at sufficiently large Q^. Indeed, the F2 - slope starts to 
deviate from the DGLAP predictions around ~ 5 — SGeV"^. 

Such a large value of this separation parameter between short and long distance processes 
leads to a natural question: "Could the experimental behaviour of the F2 - slope as a function 
of be a manifestation of the shadowing corrections or, better to say, of the saturation of the 
gluon density in DIS ?' 



The answer [|8[ g9[ gT[ is: "Yes" 



F2 - slope and SC: In Refs. I^TJ the influence of the SC on the behaviour was 

calculated using Eq.( O ) and similar formula in quark-antiquark sector. The result is given in 



Fig. 11 where the F2 - slope in ALLM parameterization |Q is also plotted. Comparison with 
the experimental data as well as with the ALLM predictions shows that the SC can give a 
plausible explanation of the effect. 

We would like to draw your attention to the fact that at fixed value of the SC do not 
lead to qualitative change of x - dependence of the F2-slope ( see lines with fixed in Fig. 11 
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Figure 10: The ZEUS data on the F2 - slope ( Caldwell plot ). 

). However, the saturation of the gluon density gives F2 oc at < Qo(^)- Therefore, 

the F2 - slope shows a decrease at small Q^. The calculations support the idea that the exper- 
immental selection of the data on the F2 - slope shows predicted Q^-dependence. However, it is 
important to stress that only data on F2 - slope at fixed energy ( or a; ) would clarify how much 
of the experimentally observed effect is related to the - dependence. It is worthwhile men- 
tioning that the ALLM'97 parameterization which describes phenomenologically the matching 
between "hard" and "soft" interactions in DIS, predicts the same behaviour of the F2 - slope 
as the SC calculations ( see Fig. 11 ). 
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Figure 11: The F2 - slope in QCD calculations taking into account SC ( GLMN ) ^1\J and 
in ALLM'97 parameterization ]5^/ , which describes the matching between "soft " and "hard" 
interactions in DIS. 
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3.3 Why these two facts are still only the indications of SC contri- 
butions? 



I hope that I have convinced you that the SC provide a natural explanation of both exper- 
imental effects. Even more, the new experimental data give a possible way to understand and 
resolve the HERA puzzle, that has been formulated at the beginning of this section. We have 
expected such effects for a long time, and the SC recieved strong support from the experiment. 

However, we cannot not make a final conclusion because both the experimental sets of data 
on a: p- dependence of the diffractive cross section, and on - dependence of the F2 - slope 
can have an alternative explanation. For example, the MRST parameterization |^ gives the 
explanation for the F2 - slope behaviour assuming the strange behaviour of the gluon structure 
function at the initial virtuality = - decreasing at low x. 

Nevertheless, we would like to stress that the SC give a good agreement with the ALLM'97 
parametrization in the wide range of and x. This provides strong support for SC which, we 
hope, will be confirmed by the future experimental data. 

We need to introduce SC in the Monte Carlo codes to be prepared for the future experiments. 



4 Matching of "soft" and "hard" photon - proton inter- 
actions 

Gribov's approach: The rich and high precision data on deep inelastic ep scattering at 
HERA , covering both low and high regions, lead to a theoretical problem of match- 

ing the non-perturbative ("soft") and perturbative ("hard ") QCD domains. This challenging 
problem has been under close investigation over the past two decades, starting from the pio- 
neering paper of Gribov ( see also ). 

Gribov suggested two stages of the 7*p interaction in any QCD description (see Fig.l2-a): 

1. The 7* converts into a hadron system (quark-antiquark pair to lowest order) well before 
the interaction with the target; 

2. The quark-antiquark pair (or hadron system) then interacts with the target. 

These two stages are expressed explicitly in the double dispersion relation suggested in Ref. 
m ( see Fig.l2-b ): 

where M and M' are the invariant masses of the incoming and outgoing quark-antiquark pairs, 
(t(M^, M'^, s) is the cross section of a qq interaction with the target, and the vertices r^(M^) 
and r2(M'2) are given by V^^W^) = R{M'^), where R^W^) is the ratio: 

R{M-) = -^-^^-^^) (26) 
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Figure 12-a Figure 12-b 

Figure 12: Gribov's approach to photon - hadron interaction: (a) the time structure of this 
interaction and (b) the double mass dispertion relation. 

which has beem measured experimentally. 

The key problem in all approaches utilizing Eq.( ^ ) is the description of the cross section 
a{M\M'^,s). 



1. We introduce Mq ^ 1 GeV in the integrals over M and M', which play the role of a 
separation parameter. For M, M' > Mq the quark - antiquark pair are produced at short 
distances { r± oc < j^^ ), while for M, M' < Mq the distance between quark and 
antiquark is too long {r± ex. > ■^), and we cannot treat this qq - pair in pQCD. 
Actually, we cannot even describe the produced hadron state as a gg - pair; 



2. For M, M < Mq we use the Additive Quark Model in which 



a{M\M'^,s) 



a^^\M\s)5{M' - M''') 
{cTtotiqN) + atotiqN)) S{M^ - M'^) 



(27) 



3. For M, M' > Mq we consider the system with mass M and/or M' as a short distance 
quark - antiquark pair, and describe its interaction with the target in pQCD. The exact 
formulae for a'^°'^'^{M'^ , M'^, s) is given in Ref. |5^), but the key property of these formulae 
that this interaction can be expressed through the gluon structure function, and it is not 
diagonal with respect to the masses, contrary to the "soft" interaction of a hadron system 
with small mass. 



Descriptions of the experimental data: Starting from the paper of Badelek and 
Kwiecinski Gribov's ideas have been implemented to describe the rich and precise ex- 
perimental data on 7*p interaction in the wide region of and x H3-It has been 



shown that such an approach is able to provide a successful description of the experimental 
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data on photon- proton interaction, over a wide range of photon virtuahties < < 100 GeV'^, 
and energies y/s (= W) < 300 GeV. The key assumption on which this aproach is based, is 
that the non-perturbative and the perturbative QCD conributions in the Gribov formula can 
be separated by the parameter Mq. The successful reproduction of the experimental data 
(see Figs. 13 ) shows that this assumption appears to be valid. It lends futher credence to 
using the additive quark model (AQM) to describe the non-perturbative contribution. 
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Figure 13-b 



Figure 13: The description of experimental data in Gribov approach, given in Ref. pl^. 

A success in phenomenological description gives us hope, that the approach based on Gri- 
bov's ideas can be used for matching of "soft" and "hard" processes for more general observables 
than the total cross section and, finally, will lead to selfconsistent Monte Carlo code for all values 
of the photon virtualities. 

Golec - Bierat and Wiisthoff approach: The natural generalization of discussed at- 
tempts to describe the matching between long and short distances physics is to use the SC 
approach. The full description based on formulae of Eq. ^ ) - type has not been developed 
yet. However, Golec - Bierat and Wiisthoff |Q suggested a simple phenomenological ap- 
proach which incorporates the main qualitative and even quantitative feathures of the general 
SC calculations. 



They use for the total cross section of photon - proton interaction the following formula 



a(7» 
where 

a{x, r±) = o-Q { 1 
R^{x) is a new scale for the DIS, namely, R^' 



(28) 



phenomenological parameterization for 



{X) 

'x) 



- e } . (29) 

\/Ql{x). Golec-Bierat and Wiisthoff use the 
Ql {x/xq)^ instead of calculation that lead 



to Fig.6. With (To = 23.03 m6,A = 0.288 and xq = 3.04 10"^. Golec-Bierat and Wiisthoff 
successfully described both the total cross sections of the photon - proton interaction as 
well as the cross sections of the diffractive production in DIS Fig. 14 show what kind of 
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description can be reached in such an approach, as well as the value of the typical momentum 
Qo{x) and effectivew Pomeron intercept in their model. 
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Figure 14: Golec-Bierat and Wiisthojf approach: ( a ) the critical line = Qq{x); ( b ) the 
effective Pomeron intercept and (c) the description of the experimental data 

This approach demonstrates that SC are a possible way of describing the matching of the 
"soft" and "hard" kinematic regions. On the other hand it gives a practical way to write the 
matching of pQCD and np QCD contribution into Monte Carlon codes. 
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Figure 15: The lego - plot of the typical LRG process of "hard" dijets production with large 
rapidity gap between them where no hadrons are produced 

5 Survival probability of Large Rapidity Gaps ( LRG) 



Definitions and general discussion: A LRG process pT] []62 | is a multiparticle production 
process in which no particles are produced in a large window in rapidity. The typical process 
of this kind is the production of two jets with high transverse momenta {pu and p2t ) and large 
rapidity difference { yi — y2 ^ 1 ) between them in which no hadrons are produced 
Fig.15 ). 



see 



hadrons [Y — yi 



p + p 

+ jet[yi,pit + LRG[yi 



(30) 



Z/2] + jet[y2,P2t] + hadron[y2-0 



Bjorken advocated that LRG processes is a unique way to measure high energy asymptotic 
at short distances. In our slang, we call this asymptotic the exchange of "hard" Pomeron, 
expressing our hope that we can calculate it in pQCD. The following observable was suggested 
5T| |]62| to be measured experimentally ( see Fig. 16 ): 

a{ dijet production with LRG ) 



gap 



(Tindusive{ dijct production 



(31) 



Unfortunately, the experimental obsevable does not measure the contribution of the 



"hard" Pomeron but has a factor |^ < \S'^\ > which we call the survival probability of LRG 



and which is the subject of this section. We will try to answer the following questions: 



1. What is the survival probability of LRG processes ? 
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Figure 16: General expression for fgap. 



2. What is the value of survival probability ? 



3. How well can we estimate this value ? 



4. What is the energy dependence of survival probability ? 

Data: We start to discuss the problem of the survival probability citing all experimental 



data on this subject given by DO|Q amd CDF||64 



DOH CDPg 

fgapiVs = ISOOGe^) 0.54 ± 0.17% 1.13 ± 0.16% 

1.85 ± 0.38 % 2.7 ± 0.9 % 

3.4 ± 1.2 2.4± 0.9 

2.2 ±0.8 ??? 



fgapiv^ = 6A0GeV] 



T3 ^ /3ap(v^ = 630 Gey ) 
sap fgap ( v/i = 1800 Gev ) 



_ (|5'P)(Vi = 630 Gey) 



(|5|2>{ = 1800 Gey) 



We can learn two important properties of (| 5* p): 



1. The value of (| 5* is small ( (| 5* < 10 % ) because all estimates show that the 
ratio of "hard" Pomeron to inclusive dijet production is about 10 - 20% ; 
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Figure 17-a 
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Figure 17: Examples of processes that lead to (| Sbremsstrahiung P)- ( Fig. 17-a ) and to (| 

SspectataA P) ( FlQ. 17-b ). 

2. (I p) decreases rapidly with energy. At least we have to blame the survival probability 
for experimentally observed fast decreasing of fgap, since the ratio of "hard" Pomeron to 
inclusive dijet production can only increase in pQCD. 



It should be stressed that only the first data on fgap in DIS |^ have appeared, this does 
not allow us to make any conclusions on the energy behaviour of (| S p) but tell us that the 
value of the survival probability, is much higher ( (| 5* p) ~ 50 — 70% ) than in hadron - 
hadron collisions. 

Two source of survival probability: It turns out that the survival probability (| 5* p) 
can be written as a product of two factor. Each of them has a different physical meaning. 

(I S H = (I Sbremsstrahlungi^y = IVl —^21) H X (| S spectator s{s) (32) 

The first factor in Eq.( ^ ) ( (| Sbremsstrahiung P) ) dcscribcs the probability that the LRG 
will not be filled by emission of bremsstrahlung gluons from partons, taking part in the "hard" 
interaction (see Fig. 17-a ). It depends mostly on the value of the rapidity gap ||6^. The 



second factor( (| Sgpectator P) ) appears to take into account the probability that no parton with 
Xi > Xi will have an inelastic interaction with any parton with x < X2 (see Fig. 17-b ). 
From Fig. 17-b one can see that this factor depends on the total energy of the process ( s ) 



p9| |]70| . Different physics behind these two factors lead to a different theoretical status for 
their calculations. For calculations of the value of (| Sbremsstrahiung P) the pQCD thechnique 
could be and has been developed Only a part of (| Sspectator P) can be controled by pQCD, 
while the estimates for most of this factor could only be made in npQCD. In pratice, it means 
that we must rely on high energy phenomenology in our attempts to provide an estimate for 

( I S spectator \ ) • 

Before discussing what we have learned about (| Sspectator P) we would like to remind the 
reader that all corrections for interactions of two groups of parton ( see Fig. 17-b ) cancel due 
to Abramovski-Gribov-Kancheli cutting rules 
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Central Production 
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Y = ln(S/)i^) 



Diffractive Dissociation 



' ► y 2 

Y = ln(S/)i^) 



Figure 18: The rapidity structure of the typical event in the Eikonal model. 



Survival Probability in the Eikonal Model: In order to understand the main features 
of the survival probability, let us consider it in the simplest phenomenological model for the 
"soft" interaction - the Eikonal model. 

• Main assumptions: 

• Only the fastest parton interacts with • ^^^f^ ^ ^ 
"wee" partons 

• Hadrons are correct degrees of freedom 

at high energy . Oversimplified final state 



Simple form of unitarity : 

2Imaei{s,b) = \ aei{s,b) |^ + Gin{s,b) . (33) 

Solution: 

n(s.b) 

aei = i[l - e 2 J ; (34) 
G,n{s,b) = l-e-^(^'''); (35) 

where Q{s, b) is an orbitrary real function. 

Simple parameterization : 

n(s,b) = —^(-)^''e~^ = iy(s) e~T^ (36) 
7rit^(s) So 
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• Analytical formulae for "soft" observable : 



atot = 2Jd%{l - e-^(^''')/2) = 27r/?2(s) [/n(^) + C- Et{~!^) ] ; (37) 

= J £h{l - e-^(^''')) = TrR\s) [ ln{v{s)) + C- Ei{-v{s)) ] ; (38) 

<yei = cTtot - cT,n = 7^R\s) [/n(^) + C + Ei{-v{s)) - 2Et{-^-^) ] . (39) 

A simple formula for survival probability : 

< I ^spectator- I >- Jcmrnib) ' ^ ^ 



with r^,(6) = 

Analitycal form of <\ S spectator P> ■ 



Q |2^_ a{sh[a{s),iy{s)] 

^1 '^spectator \ ^ r / v / 



\V 



with a{s) = 



Survival Probability in the Eikonal Model ( comparison with experiment ): 

To get more reliable estimates for the value and energy dependance of <| S spectator P> we 
tried to extract all parameters directly from the experimental data [|70|. The first observation 
( see Eqs. ( |3^ ) - ( ^ ) ) is that the ratio 

Relil^) = — (42) 
C^tot 

is a function of the only parameter u. Therefore, we can use the experimental data on Rei to 
fix u at particular value of energy ( see Fig. 19 ). 

To fix the value of a{s) we used the data on J/\l' production |Tl|] which shows two differents 



slopes in t-dependence for elastic and inelastic processes ( see Fig. 20-a ). These data can be 
used to extract the value of Rjj = SGeV'"^ ( see Ref. ). 

Using "soft" phenomenology for -R^(s) we obtained prediction for the survival proba- 
bility given in Fig. 20. One can see that the Eikonal model reproduces both the experimental 
value and the energy dependence of the survival probability. 

More general approach to calculation of (| Sgpectator |^)- 

In Ref. [0 a more general approach has been developed to calculation of <| Sgpectator P> 
which based on the following assumptions: 



• Only the fastest parton interacts with the target ; 



• Hadrons are not correct degrees of freedom ; 
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Figure 19-b 



Figure 19: Experimental data ( Fig. 19-a ) and the Eikonal model predictions ( Fig. 19-b ) for 

Rel = ^ ■ 




Figure 20-a Figure 20-b 

Figure 20: J/Psi production (Fig. 20-a ) and the survival probability <\ Sgpectator P> the 
Eikonal model ( Fig. 20-b ) 
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^ (7dd{M>Mo) ^ . 

• Oversimplified final state: central production ( uniform in rapidity) and elastic + diffrac- 
tive production with small masses. 

One can see that all the above assumptions are very close to the Eikonal model, and the 
difference actually is only in the way of how the diffractive production processes are treated. 
In the model of Ref. |^ we assume that the correct degrees of freedom at high energy are 



not hadrons but some different states which are described by wave functions and \l/2. Our 
interaction matrix is diagonal 

< X ^'„|T|^^„/ ^-^m'> = ^n,m(s, h) 5{n - nl) 5{m - m') , (43) 

and only for An,m ( not for hadrons ) unitarity has a form: 

2ImAn,m{sA) = \Alm{sM)? +Gl^{sM). (44) 
Produced hadrom and a diffractive state have the following wave functions: 

"^hadron = « ^1 + /5 ^2 ; (45) 

■q/j^ = -/? ^1 + a $2 • (46) 

Therefore, three channels were taken into account in this model, namely, elastic scattering, 
single and double diffractice production while only elastic rescatterings were included in the 
Eikonal model. 

Using the Eikonal - like parameterization ( see Fig. 21 ), the three channel model leads to 
the prediction for the value and energy dependence of (| Sspectator P) given in Fig. 22. One can 
see that this model can reproduce the experimental data including the energy dependence. 

Summarizing the experience dealing with three channel model, we can conclude that: 



• Accuracy of the three channel model, in principle, is much better than in the Eikonal 
Model, but still we do not know what to do with diffraction in the region of large mass ( 
M ^ s ) ; 

• The scale of SC is not given by Rd = ^ ' rather by the separate ratios 

^' (Ttot (Ttot (Ttot 

• The small value of the survival probability as well as its strong energy dependence appear 
naturally in this approach ; 

• The parameters that have been used are in agreement with the more detailed fit of the 
experimental data ; 



Questions to experimentalists: 
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Figure 21: T/ie Eikonal - like parameterization in the three channel model. 

1. What is Rei at ^/s equal to 630 GeV ? 1800 GeV ? Tevatron runll energy ? 

2. What is the value for Rsd at ^/s equal to 630 GeV ? 1800 GeV ? Tevatron runll 
energy ? 

3. What is the value for Rdd at equal to 630 GeV ? 1800GeV ? Tevatron runll 
energy ? 

4. What is the value of the single diffraction cross section in the region of large masses 
at y/s equal to 630 GeV ? 1800 GeV ? Tevatron runll energy ? 

Q & A: 

Q: Have we developed a theory for (| Sgpectator P) ? 
A: No, there are only models on the market . 

Q: Can we give a reliable estimates for the value of (| Sgpectator P) ? 

A: No, we have only rough estimates based on the Eikonal - type models. 
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Figure 22: The value ( Fig.22-a ) of survival probability <| Sgpectator P> o^t yfs = W = 1800 Gel/ 
and the ratio <\ S{W = GAOGeV) \^> / <\ S{W = 1800 Gel^) p ( Fzg.22-b ) versus p in 
the three channel model. 



Q: Can we give a reliable estimates for the energy behaviour of (| S spectator ? 

A: No, but we understand that (| Sspectator P) could decreases steeply with energy. 

Q: Why are you talking about (| S spectator P) if you can do nothing ? 

A: Because: 



• Dealing with models we have learned what questions we should ask experi- 
mentalists to improve our calculations; 

• We learned what problems we need to solve theoretically to provide reliable 
estimates; 

• We understood what kind of questions could be answered in LRG experiments; 

• We are on the way to an understanding in what experiments our accuracy is 
enough, to get interesting information on the high energy scattering amplitude 
at short distances. 

Survival Probability in DIS: I hope, that I have explained to you that we do not 
have a consistent theoretical approach to the calculation of (| Sspectator P) for hadron - hadron 
collisions. The natural question to ask is : can the theoretical situation be better for DIS where 
the "short" distances can give the major contribution? 
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Figure 23: The three channel model predictions for double diffractive dissociation ( Rdd 
(^dd/ O'tot- 
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Figure 24-a 



Figure 24-b 



Figure 24: Lego - plots for DD ( Fig. 24-a) and one-side dijets production in DD ( Fig. 24-b ) 
%n DIS. 



To answer this question let us consider first two examples of the LRG processes in DIS: the 
diffractive dissociation ( DD ) and one-side dijets production in DD ( see Fig.24 ). 

Fig. 25 shows the parton interaction which can occur for one-side dijets production in DD. 

Using the main properties of the parton cascade ( mostly the AGK - cutting rules ||6^ ) we 
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Figure 25: Picture of parton interactions for one-side dijets production in DD in DIS. 



can make three observation on the survival probabihty for this process: 



1. Survival probability for one - side dijets depend on and x^, but does not depend on 
Pt and Xi ; 

2. Survival probability is not big in comparison with LRG survival probability in proton - 
proton collisions ; 

3. For large Q"^ the survival probability can be calculated in pQCD . 



In the Glauber - Mueller approach the cross section of virtual photon - proton interaction 
can be written in the form []73| []73|: 

<1*P) = j j d\dz \^iQ';z,r,)\'aix,r',) ; (47) 

with 

^{x,r^) = 2 J d%, [ 1 - e- "-^"^^ ] ; (48) 
where ^ is the wave function of the virtual photon and 

^{rlxB-M) = ^r,^^BG(xB,|)-^e-*. (49) 



i?2 ^ 2B'^^{xb) where 

e 



^(^ = 0) 

and this parameter is closely related to non pertutbative physics. At the moment we can use 
the experimental data on the slope of the DD processes in DIS. It should be stressed that 
factoring out the bt - dependence can be justified in QCD using the factorization theorem |jT5[ . 
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The formula for the survival probability is a direct generalization of Eq. ( |^ ) in the Eikonal 
model, namely: 



\\ >^ spectator \ / 2b2 



(50) 



The calculations , using Eqs. ( ^ ) - ( 0, shows that the survival probability is much 
less in DIS than in hadron - hadron collisions ( see Fig. 26, where = S"^ is plotted for the 
transverse polarized induced photon reactions with LRG ). 
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Figure 26: The survival probability (| Sspectators P) = D"^ in DIS versus x ( Fig. 26-a ) and 
( Fig. 26-b ). 



Multiparticle production in DIS: Shadowing corrections or, in other words, the rescat- 
terings of quark - antiquark pair in a target lead to a new predictions for the processes of 
multiparticle productions. Namely, it gives probabilities for production of several parton show- 
ers. For every parton shower we can use the Monte Carlo code based on the evolution equations, 
while the probability for configuration with n - parton showers can be calculated, using the 
following simple formula: 

a„ = y d\dzd%t\m{Q\rt,z)\^ —e-"" (51) 

Fig. 26 shows the first three terms in our decomposition of the multiparticle production in the 
parton showers: diffraction dissociation, one parton shower and two parton showers productions. 
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Figure 27: The first three terms in parton showers decomposition. 

For better understanding, in Fig. 26 the probability of each configuration in the hmit of small 
Q is also given. 

I am firmly believe that such a decomposition will be very useful for writing Monte Carlo 
programs. 



6 Conclusions 

In this talk I tried to give you a picture of recent progress in low x physics. I chose the subject 
which, I hope, would be useful for Monte Carlo experts. Let me summarize in short the present 
situation in four subjects that I have discussed here: 

1. We are on the right track in understanding of the next-to-leading BFKL Pomeron, but, 
at the moment, we cannot guarantee any calculations of the parameters of the BFKL 
Pomeron; 

2. First experimental indications of the strong shadowing corrections have appeared in en- 
ergy behaviour of the diffraction dissociation cross sections in DIS and in the - be- 
haviour of the F2 - slope ( Caldwell plot ). We hope, that these data will stimulate the 
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new experimental systematic search of SC, as well as the creation of the new Monte Carlo 
codes, that will include the SC and physics, induced by them; 



3. Considerable progress, based on Gribov's ideas, has been achieved in describing the match- 
ing between the non-perturbative "soft" processes, and the perturbative "hard" ones in 
photon - proton interactions. We firmly believe that this progress will be useful both for 
future non-perturbative approaches in QCD and for creating new Monte Carlo program, 
which take this matching into account; 

4. We can estimate the survival probability in hadron - hadron collisions only in rather 
primitive models, which can reproduce the basic experimental data on it. The prediction 
for the survival probability in DIS are more solid from theoretical point of view, but 
the experimental information is so poor that we cannot compare the main features of 
our calculations with the data. We think that time has come to include in Monte Carlo 
codes the decomposition of multiparticle production process in the series of multi parton 
showers production, which is closely related the the value of the survival probability. 
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